Persistence With Partial Survival 
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We introduce a parameter p called partial survival in the persistence of stochastic processes and 
show that for smooth processes the persistence exponent 9{p) changes continuously with p, 9{Q) 
being the usual persistence exponent. We compute 6(p) exactly for a one-dimensional deterministic 
coarsening model, and approximately for the diffusion equation. Finally we develop an exact, 
systematic series expansion for 6(p), in powers of e = 1 — p, for a general Gaussian process with 
finite density of zero crossings. 



Recently considerable theoretical and experimental ef- 
fort has been devoted to understanding first-passage 
statistics in nonequilibrium systems. These include the 
Ising, Potts and time-dependent deterministic Ginzburg- 
Landau (TDGL) models undergoing zero-temperature 
phase-ordering dynamics [0-^, the diffusion equation 
with random initial conditions the global magneti- 
zation undergoing critical dynamics , several reaction- 
diffusion systems Q , fluctuating interfaces Q and a ran- 
domly accelerated particle Typically one is inter- 
ested in 'persistence', i.e. the probability P[){t) that, at a 
fixed point in space, the stochastic process (such an Ising 
spin or the diffusion field) does not change sign up to 
time t. In the examples mentioned above, this probabil- 
ity decays as a power for large time t, Po{t) ~ , where 
the persistence exponent 9 is nontrivial due to the non- 
Markovian nature of the process in time at a fixed point 
in space. This exponent has recently been measured ex- 
perimentally in a 2-d liquid crystal system and also 
for 2-d soap froth [|l2| and breath figures |l3|] . The the- 
oretical computation of 9 however, despite a few exact 
and approximate results, remains a major challenge. 

Even for the simple diffusion equation, dtcj) = 
starting from random initial configuration, the exponent 
9 is known only numerically and within an independent 
interval approximation (IIA) though there is a re- 
cent conjecture for an exact 9 that remains to be 
proved. The IIA result, though in excellent agreement 
with numerical simulations, is hard to improve system- 
atically. The central result of this Letter is to derive a 
systematic series expansion for 9 in terms of a suitable 
expansion parameter. This expansion is exact order by 
order and when truncated at second order already gives 
good results for the diffusion equation. But this exact se- 
ries expansion technique is more general and goes beyond 
the diffusion equation. We show that it can be applied to 
compute the persistence exponent, order by order, for a 
wide class of stochastic processes which includes the dif- 
fusion equation, random acceleration and the 1-d TDGL 
model as special cases. 

Our result is also useful for a related problem which 
has wide applications in diverse fields ranging from 



information theory to stock markets and oceanogra- 
phy. Consider a stochastic Gaussian stationary process 
X{T) characterized completely by its two-time correla- 
tor, (X(O)X(T)) = /(T). The process X{T) can be 
used to model, e.g., the current in an electrical circuit 
or the price of a stock. Given /(T), what is the prob- 
ability Po{T) that the signal stays above (or below) a 
certain level , say zero, up to time T? This problem has 
been studied for many years p6| , p^ and it is known that 
ff |/(T)| < l/T for large T, then Po{T) - exp(-6ir) 
for large T |16|. The exponent 9 depends quite sensi- 
tively on the full function f{T) and is very hard to com- 
pute for general f{T) Q. For a Markov process, where 
/(T) = exp(-AT) for aU T, it is known that 9 ^ X @. 
However for non-Markov processes, where f{T) is not a 
pure exponential, very little is known. Only recently a 
perturbation theory result for 9 was developed for pro- 
cesses close to Markovian [§,0 . 

The persistence problem for the diffusion equation in 
d-dimensions can be exactly mapped to a Gaussian sta- 
tionary process, with f{T) — [sech( T/2)]'^/^, b y identi- 
fying T = Int and X(T) = (j){x,t) / ^ {4>^\x,t)) §. The 
probability of no zero crossing then decays as Po{T) ~ 
exp(— 6'T) — t^^ . The series expansion technique that 
we develop below can be used to compute 9 for arbitrary 
/(T) as long as /(T) ~ 1 - aT^ -h ... for small T. Such 
Gaussian processes are called smooth as they have a finite 
density of zero crossings, p = /"(0)/7r 

The key strategy underlying our technique is to first 
generalize the usual persistence problem by introducing a 
partial survival factor p as follows. The usual persistence, 
say in the diffusion equation, is the fraction of points in 
space where the diffusion field has not changed sign even 
once up to time t. One way to compute this is to start 
with a random initial configuration of the field and put 
a particle at each point in space to act as a counter. At 
subsequent times, whenever the field changes sign at any 
point, the particle there dies. The persistence is simply 
the fraction of particles still surviving at time t. We now 
generalize this by assigning the rule that whenever the 
field changes sign at a point, the particle there survives 
with probability p and dies with probability (1 — p). We 
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then compute the fraction of particles, P{p,t) left after 
time t. Thus, p = corresponds to usual persistence 
P{0,t). A somewhat similar generalization was recently 
studied in the context of "adaptive persistence" problems 

This generalization has several implications. It is easy 
to see that if P„ (t) denotes the probability of n zero cross- 
ings in time t of the underlying single site process, then 
P{p, t) is simply the generating function, 

oo 

P(p,i) = ^p"F„(t). (1) 

ri=0 

For p — P{0,t) ~ Po{t), the usual persistence, decay- 
ing for large t as t~^^^\ In the other limit, p — 1, the par- 
ticles always survive: P{l,t) = 1, implying 9{1) — 0. It is 
interesting to analytically continue Eq. to negative p. 
In fact, for p = —1 this is simply the autocorrelation 
function, P(-l,t) = A{t) = (sgn(0(x, 0)) sgn((/)(x, t)), 
which decays as t~^/'^, where A is a well-studied expo- 
nent in phase-ordering systems ||2^. In fact, we show 
below that for smooth processes, P{p,t) ~ t-^(p) for 
large t where the exponent d{p) depends continuously 
on p as p varies from —1 to -1-1. Moreover, the quantity 
Apit) — P{—p, t) I P{p, t) is just the autocorrelation func- 
tion averaged only over points with surviving particles, 
when the survival probability is p. So if Ap{t) ~ 
we have Ap — 0{^p) — 0{p)- this generalization thus puts 
both the autocorrelation and the persistence exponents 
as members of a wider family of exponents. 

We first establish the continuous dependence of 9{p) 
on p for smooth processes by computing 6{p) exactly 
for a non-Gaussian process, namely the 1-d determin- 
istic TDGL model, and then approximately within IIA 
for the diffusion equation. We then proceed to compute 
9{p) for any smooth Gaussian stationary process by ex- 
panding around p = 1. This series expansion result for 
6{p) in powers of e = 1 — p is exact order by order. 

If a system, such as the Ising model, is quenched from 
a high-temperature disordered phase to zero tempera- 
ture, domains of 'up' and 'down' phases form and grow 
with time. The evolution of the order-parameter field 
(f) can be modelled by the deterministic TDGL equa- 
tion, dt4) — — V'{4)), where V{(t)) is a symmetric 
double well potential with minima at = ±1. In 1-d, 
at late times the system breaks up into alternate 'up' 
and 'down' domains and coarsens by successively elim- 
inating the boundaries of the smallest domain, i.e., by 
flipping the signs of <j) simultaneously at all points in- 
side the smallest domain The density of persistent, 
or 'dry' parts where (j) has not changed sign then scales 
as ~ (/)~^^°-' where {I) is the average length of growing 
domains, which serves as 'time' in this problem. The 
exponent 0(0) was computed exactly by noting that the 
dynamics does not generate correlations between neigh- 
boring domains 0. We now introduce the partial sur- 



vival factor p in this dynamics. 

We start with a random distribution of intervals or 
domains and assign a particle to each point in space. 
The dynamics merges the smallest interval Imin with its 
two neighbors Ii and I2 to make one single interval /. 
The lengths and the 'dry' part d{I) (i.e, the num- 
ber of live particles in the interval /) evolve as, 1{I) — 
l{h) + Kl2) + l{Irmn) and d{I) = d{h) + d{l2)+pd{Imin)- 
Thus the only difference from the calculation in Ref. |^ 
is the p-dependent term in the dry part. The rest of the 
calculation is similar to that in Ref. |^ and we just out- 
line the method without details. One writes down the 
evolution equations for the number of intervals of length 
I and the average dry part carried by such an interval, 
and one solves exactly for the associated scaling func- 
tions by taking Laplace transforms. Demanding that the 
first moments of these scaling functions are finite gives a 
transcendental equation for 0{p), 




dte-' t-^-^ [(1 -p){l-t- e~*)e'-(*) 



+26{l +p)t + e(l- p)i2e-'-W] = 0, (2) 

where r{t) = —l~^'^=i{~tY'' l'nn\, 7 being Euler's con- 
stant. Clearly, for p = 1 one gets 9{1) = from the 
above equation as expected. For p = 0, it reduces to the 
equation for 9{0) as obtained in Ref. j^. Forp = —1, one 
recovers the equation for 9{—l) = A of Ref. Fig. 1 

shows 9{p) as a function of p for — 1 < p < 1, obtained 
by numerically solving Eq. (0). 
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FIG. 1. Dashed line: The exponent e{p) for the 1-d TDGL 
model, obtained from Eq. (^. Solid lines: The IIA estimates 
for 6[p) for the diffusion equation in (bottom to top) 1, 2 and 
3 dimensions. 

We now turn to the diffusion equation, dt4> — V^0, 
starting from a random initial configuration. We first 
carry out a numerical simulation to compute P{p,t) for 
finite p following the procedures of Ref. ||^. Fig. 2 
shows the asymptotic decay of P{p,t) with i on a log- 
log plot for p = and p = 0.5 in 1-d. Clearly the 
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exponents are quite different. For example, for p = 0, 
e{0) = 0.1207 ± 0.0005 as in Ref. g but for p = 0.5, we 
find 9{l/2) = 0.0588 ± 0.0005. Unfortunately we have 
not been able to compute 9{p) exactly. However, the IIA 
used in Ref. ||^ to compute 9(0) can be easily extended 
to compute 9{p) very accurately for all p as we now show. 
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FIG. 2. Time-dependence of the generalized persistence 
probability, P{p,t), from numerical simulations. Circles: the 
1-d diffusion equation for p = (filled) and p = 0.5 (open). 
Squares: the 1-d Glauber model at zero temperature for p = 
(filled) and p = 0.5 (open). The lines are guides to the eye. 

Consider the normalized process X(T) ~ 4>/{4>'^) as 
a function of T = \at. The zero-crossing events of 
are the same as those of X, but X{T) is a stationary 
Gaussian process characterized completely by its two- 
time correlator, /(T) = (X(O)X(T)) = [sech(T/2)]'^/2. 
With a non-zero survival factor p, the fraction of live 
particles after time T is then, P(p,T) = Y.'^^qP^ Pn[T) 
and will decay at late times as exp (— 6'(p)T) — 
To evaluate Pn{T), the probability of n zero crossings 
by X in time T, we note that the Laplace transforms, 
P„(s) = dTeyii>{-sT)P„{T), were evaluated in Ref. 

using IIA, i.e., assuming that successive intervals be- 
tween zero crossings of X are statistically independent. 
Using these results from Ref. §, and carrying out the 
sum over n, gives P{p,s). Since P{p,T) ~ exp(— ^?(p)r) 
for large T, P{p, s) will have a simple pole at s = —9{p). 
Using this, we finally get 9{p) as a solution of the equa- 
tion, 

= eiT^I^, ( 1 + - r dT exp(0T) 

1 + P \ d I TT Jo 

X sin-i[sech'^/2(rp/2)]}. (3) 

The solution is plotted in Fig. 1 for d = 1,2,3. We note 
in the two extreme limits, p = 1 and p = — 1, the IIA 
gives 9{1) = and 9{—l) = d/A, which are exact. For 
intermediate values of p, the IIA results are in excellent 
agreement with numerical simulations. For example, for 



p = 1/2 the IIA gives 9iia = 0.05823044. . ., compared 
to 9s^rn = 0.0588 ± .0005. 

Having established the continuous p-dependence of 
9{p) for two smooth processes, we now derive an exact 
series expansion of 9{p) near p = 1 for a general smooth, 
Gaussian, stationary process X{T), characterized by its 
two-time correlator f{T). The basic idea is straightfor- 
ward. We start with the definition (|l|) of P(p, t) as a gen- 
erating function. Writing p" — cxp(nlnp) and expanding 
the exponential, we obtain an expansion in terms of the 
moments of n, the number of zero crossings: 

lnP(p,T) = f;ii^(n'-)„ (4) 

r=l 

where {n^}^ are the cumulants of the moments. Using 
p = 1 — e, we express the right-hand side as a series in 
powers of e. Since P(p, T) is expected to decay for large 
T as exp(— 0(p)T), we obtain a series expansion of 9{p) 
by taking the limit, 

oo 

9{p) = - hm - In P(p, T)^Y. (5) 

1 — ^00 !_ 

r=l 

The coefficients a^'s involve the cumulants. 

Fortunately the computation of the moments of n is 
relatively straightforward, though tedious for higher mo- 
ments. For example, the first moment (n), i.e. the ex- 
pected number of zero crossings in time T, was com- 
puted by Rice |l|]: (n) = T^--/"(0)/7r, implying ai = 
\/~f"{0)/n. The second moment, (n^), was computed 
by Bendat ||2^. Using this result and after some algebra 
we have computed the coefficient 02 , which already looks 
complicated. We just quote the final result here (details 
will be published elsewhere ]2^]): 

a2^\ I [S{oo) - S{T)]dT, (6) 

where S{T) is given by 



with H = M2i/^/M^2 - ^24- The Mij's are the cofac- 
tors of the 4x4 symmetric correlation matrix C between 
4 Gaussian variables [X{Q), X{0), X(T), X{T)]. The el- 
ements of C can easily be computed from the correlator 
/(T). For example, Cn = (X(O)X(O)) = /(O), Cm ^ 
{X{Q)X{T)) = /'(T), C24 - {X{Q)X{T)) = -/"(T) and 
so on. 

Although these expressions look complicated, in many 
cases the function S{T) can be evaluated explicitly 
and the integral for 02 can be performed analytically. 
For example, for 2-d diffusion equation, where /(T) = 
sech(T/2), we get 
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Keeping terms up to second order and putting e = 1 
(in the same spirit as e expansion in critical phenom- 
ena) gives 9(0) = (7r + 4)/47r2 = 0.180899..., just 3.5% 
below the simulation value, Osim = 0.1875 ± 0.0010 [||. 
Note that though the IIA estimate, 9ija = 0.1862 is even 
closer to the simulation, it can not be improved systemat- 
ically. The series expansion estimate, on the other hand, 
can be improved systematically order by order. 

We have also computed, for the first time, the third 
moment (n^) for a general smooth correlator f{T). We 
then use this to compute 03 . The expressions involve the 
elements of a 6 x 6 correlation matrix and are not partic- 
ularly illuminating ||2^, so we skip the details here. As 
an example, we computed the series up to third order 
for the random acceleration process, <Px/dt^ = rj irj is 
a Gaussian white noise) which can be transformed to a 
Gaussian stationary process with f{T) ~ [3 exp(— T/2) — 
exp(-3r/2)]/2 §. We find 



(9) 



Putting e = 1, we get to third order, 0(0) = 
0.271835775... which should be compared to its exact 
value 0.25 [|l^. We note that the series oscillates around 
the exact value 0.25 as the order increases. 

We note that the series expansion will fail for non- 
smooth Gaussian processes, whose moments of zero cross- 
ings are not finite. As an example, consider ordinary 
Brownian motion, dx/dt = rj, which can be mapped 
to a stationary Gaussian Markov process with correla- 
tor /(T) = exp(— T/2) using the change of variables 
discussed before. For this process it is well known 
that the moments of zero crossings are infinite: if 
the process crosses zero once, then it crosses again in- 
finitely many times immediately afterwards |^ . Thus 
only the n = term contributes to the sum (^), giv- 
ing P(p,T) « P(0,T) - exp(-r/2) = for large 
t. Thus e{p) = 1/2 for all < p < 1, except at p = 1 
where 9(1) = 0. Since 9{p) is discontinuous at p = 1, no 
expansion around p = 1 is possible. 

The same conclusion holds for the T — Glauber dy- 
namics of the Ising model. In this case, the usual per- 
sistence exponent 9{0) was recently computed exactly in 
1-d 1^ and approximately in higher dimensions The 
exact value in 1-d is 9{0) = 3/8 Even though the spin 
Si(t) at a given site i is no longer a Gaussian process, it 
is non-smooth nevertheless, i.e. if a spin flips once, it usu- 
ally flips many times immediately afterwards. This fact 
can be tested easily by computing the exponent 9{p) for 
nonzero p. In Fig. 2 we show the asymptotic dependence 
of P{p, t) on t on a log-log plot for p = and p = 1/2 for 
the 1-d T = Glauber model. In contrast to the diffu- 
sion case, the asymptotic slopes are the same and given 



by 0.375 ± 0.002. We have checked this fact for other 
values of p, and conclude that 9{p) is independent of p 
for < p < 1 (26|, while clearly 9(1) 0. Thus the p de- 
pendence of 9{p) provides important information about 
the nature of the smoothness of the underlying stochastic 
process. 

We thank Deepak Dhar for useful discussions. 
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